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Abstract 

This paper presents the dynamic behavior of a Free- 
Piston Stirling Engine/linear alternator (FPSE/LA) 
driving a single-phase fractional horse-power induction 
motor. The controllability and dynamic stability of the 
system are discussed by means of sensitivity effects of 
variations in system parameters, engine controller, 
operating conditions and me ch a n ical loading on the 
induction motor. The approach used expands on a com- 
bined mechani cal and thermodynamic formulation 
employed in a previous paper [Ref. 7], The application 
of state-space technique and frequency domain analysis 
enhances understanding of the dynamic interactions. 

Engine-alternator parametric sensitivity studies, 
sim ilar to those of the previous paper, are summarized. 
Detailed discussions are provided for parametric varia- 
tions which relate to the engine controller and system 
operating conditions. The results suggest that the con- 
trollability of a FPSE-based power system is enhanced 
by proper operating conditions and built-in controls. 


Introduction 

The potential for long life and reliability has culmi- 
nated in the consideration of the Free-Piston Stirling 
Engine/linear alternator (FPSE/LA) for space power 
applications. This realization has generated widespread 
interest in the modelling and analysis of the FPSE, for 
better understanding of its performance characteristics. 
Earlier papers [Refs. 1 to 6] have dealt mainly with the 
FPSE and its dynamics, with minimal treatment of, or 
no reference to the electrical portion of the subsystem, 
namely, the LA. Some of the more recent papers [Refs. 
7 t 8 ] have included a more detailed analysis of the LA 
The results in Reference [7], on a FPSE/LA connected 
to a resistive load, provide useful information on the 
engine-altemator-load interactions, and may be used to 
guide system design and/or experimental efforts. In the 
potential application of the FPSE/LA in space power 
systems, there is a high probability that there will be 
more dynamic than static loads. The introduction of a 
dynamic load into any existing power source-load 
interactions complicates the system sufficiently to 
warrant an assessment of the resulting system behavior. 


This paper describes the control and dynamic 
stability analysis of an FPSE/LA feeding into a single- 
phase, fractional horse-power induction motor which is 
driving a comparably-sized fan load. The engine/LA is 
the Space Power Research Engine, SPRE, [Ref. 9] at 
the NASA Lewis Research Center, which was used in 
Reference [7]. The analysis is based on a combination 
of dynamic and thermodynamic analysis of the Stirling 
engine [Refs. 7, 10], circuit theory of the electrical 
system and state-space technique in modern control 
theory [Ref. 11]. Dynamic stability is determined from 
computation of eigenvalues of the linearized system 
matrix, using the MATLAB software [Ref. 12]. Control- 
lability of the integrated system is evaluated by paramet- 
ric sensitivities on the eigenvalues, to show the potential 
of a FPSE/LA-dynamic load application in space power 
system. 


System Representation 

Figure 1 depicts the equivalent circuit of the system 
of study. It comprises a Free-Piston Stirling Engine 
coupled with a linear alternator (FPSE/LA), connected 
via control circuitry to a single-phase induction motor. 
A connected fan is driven by the induction motor- 
generated electromagnetic torque. The parameter Cp is 
a tuning capacitor for optimum power transfer to the 
motor-fan load. The term Rp denotes a parasitic load 
for contr olling the system through maintenance of power 
balance between the generation source and the motor- 
fan load. The resistor R c is inserted to determine its 
effectiveness in controlling the system during load-side 
fault. 

The parameters r s and Ljg represent the induction 
motor stator winding resistance and leakage inductance, 
respectively, of the motor. The magnetizing branch of 
the air-gap is denoted by L m . The parameters Lj r and 
r r are the rotor leakage inductance and winding resis- 
tance, respectively, referred to the stator side. The rotor 
slip, S, is a measure of how much the frequency of the 
rotor-induced field lags the stator field. 

Equations of the system dynamic behavior are 
derived, based on Fig. 1 and the schematic representa- 
tion of the FPSE in Fig. 2. 
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Development of System Equations 


Equations of LA-Dynamic Load Subsystem 


The key FPSE equations are summarized. The 
equations of the electrical portion of the system are 
developed next. The mechanical-thermodynamic and 
electrical equations together are then cast in a state- 
space notation. Dynamic stability analysis is then 
performed. 

Summary of FPSE Equations 

The dynamic equations of the FPSE relate to the 
single cylinder representation in Fig. 2, and are based on 
the assumptions below: 

(1) Schmidt’s thermodynamic analysis is evoked such 
that: 

(a) the displacer and piston motions and the 
working space pressure are sinusoidal. 

(b) the working fluid obeys the ideal gas law, 
and expands and compresses isothermally. 

(2) The working space gas pressure is spatially 
constant but time-variant. 

(3) The bounce space pressure balances the average 
working space pressure. Hence, average positions of the 
power piston, displacer and cylinder casing are station- 
ary. 

The detailed development of the pertinent equations 
is described in Reference [7]. The engine mechanical 
and thermodynamic properties are included. The 
resulting dynamic equations are stated in Eqn. (1), 
where the definitions of the terms are summarized in 
Table 1. 

[ m D 0 *D = _ C 11 C 12 *D . K 12 X D ,[o 1 

[0 m p Xp c 2lC22 X p K 21 K22 X p [ F EPj 


The loop equations of loops 1 through 4 in Fig. 1. 
are summarized in Eqns. (2) through (6), respectively, 
using Kirchhoffs voltage law. The notation / represents 
df 

dt ' ^ 

(R*+ Rp+ Re) ig + kma) * g ~ ^ p* L + * c g * dSCp^p 

-Rpig* fr** V'L* < L ls + \ r 0.SLj r OSLj h -0 (3) 

— 0JLq|1^ + 0.5 y+ 0.5(L m + L b ) iy -0 (4) 

- O-SlJl + 0.5 r£ i b + 0.5(L m+ L^) ? b - 0 (5) 

(2*S) 

V C T - k 1 C T (®) 

/ \ • 

The term ^ is the voltage generated by the 

piston-induced flux change. The loop currents L and ^ 
relate to the forward and backward rotating fields in the 
rotor [Ref. 13]. The electrodynamic differential equa- 
tion (7) relates to the dynamics of the induction ma- 
chine-fan load. The 

Jw r = T e - kw r 2 (7) 

parameters J and w r are the respective inertia and 
rotational velocity of the induction machine rotor, k is 
the rotational constant of proportionality for the fan. 
The nominal value of k is chosen to correspond to the 
033 hp rating of the induction motor. The terms 
and R b are the resistive components of the impedances 
Zy and Zfc ofthe forward and backward rotating circuits 
in Fig. 1. The derivation of the electromagnetic torque 
T e , developed by the induction machine yields Eqn. (8), 
where w cx is the excitation angular velocity of the stator 


where 


[ C 11 C 12 


(C do +C T ll) teop +C Tl2> 

[ c 21 C22 


(Cpp+Cj2l) (Cpp +C722), 

F K n *.2' 


{K D d+ k tii) + 

F K 21 *22 


{ K PD + KT2l) (Kpp + K722) 




T t = CONST2*i L 2 = ((1-S) (Ry - R*) Pp/4w ex ) i L 2 (8) 

winding. The slip S of the forward field relates to the 
relative motion of flux and rotor conductors, with 
subsequent induced rotor voltage of slip frequency S/ cr 
The parameter P p is the number of poles on the induc- 
tion machine. 


System Dynamic Stability 

The combined dynamic equations of the FPSE and 
the circuit and electrodynamic equations of the electri- 
cal-mechanical load subsystems may be linearized for 
dynamic stability analysis. State variable representation 
[Ref. 11] of the linearized equations enable easy compu- 
tation of the system eigenvalues. 
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State-Space Formulation 

The no nlin earities in Eqns. (7) and (8) are removed 
by linearization about the steady-state operation. 
Neglect of second order incremental terms results in 
Eqn. (9). 

jAw r = p.CONSm.ijJ Ai L -2 kw ro Aw r (9) 

The system state variable vector in terms of incremental 
changes may be defined as in the transposed (T) Eqn. 
( 10 ). 

[X) T s [Xj Xj X 3 X 4 X5 X, X7 Xg x, X 10 ] 

= IaX D AX d AXp AXp Ai g AV cT Ai L ti f Aig AW r ](10) 

Substitution of Eqn. (10) into Eqns. (1) through (6) and 
(9) results in the state-space equation (11), where 

[X] = [A][X] (11) 

the system [A] matrix is defined in Eqn. (12). The 
matrix entries in relation to the system parameters are 
summarized in Table 2. 

Discussions of the dynamic stability analysis on the 
eigenvalues of the [A] matrix follow. 

Results of Sensitivity Analysis 

Simulation of the system behavior under varying 
parameters and operating conditions is performed, using 
the MATLAB software [Ref. 12]. The results of the 
simulation appear in Figs. 3 through 10. The base or 
reference system results of Fig. 3 are obtained from the 
use of the nominal parameters in Table 1. The FPSE 
and LA subsystems are the same as those used in 
Reference [7], namely, the SPRE engine [Ref. 9]. 
Subsequent figures show the particular parameter being 
varied, while the remaining parameters are kept at the 
nominal values. The physical meaning of the eigenvalue 
location in the complex plane is given by the modal 
response, such as illustrated in Fig. 4 of Reference [7]. 

Fig. 3 shows the electrical, piston and displacer 
roots El, MP and MD, respectively, of the coupled, and 
the corresponding roots E1U, MPU and MDU of the 
uncoupled FPSE and LA plus associated control param- 
eters Cj and Rp. Similarly to the presentation in 
Reference [7], these roots yield the same information 
with regard to the design (105 Hz) and operating (101.1 
Hz) frequencies, respectively, of the SPRE engine. Also, 
Redlich and Berchowitz’s observation [Ref. 4] of the 
piston root, MPU, in the right-hand plane is confirmed. 


This root is indicative of an oscillator with amplitudes 
which tend towards its mechanical limits. Furthermore, 
the coupled system’s reduction in the margin of operat- 
ing-point stability of the electrical subsystem is evident. 

The presence of the induction machine load yields 
additional three electrical roots E2, E3 and E4. These 
roots are associated with the main field of the stator, 
and the backward and forward rotating fields of the 
rotor, respectively. The real-only nature of these roots 
signifies that, for the small horse-power (1/3 hp) ma- 
chine in question, any inherent electric transients are 
typically highly damped. In Fig. 1, the portion to the 
right of the Rp branch is the steady-state equivalent 
circuit of the angle-phase induction motor. The absence 
of machine transients from the circuit eliminates any 
conjugate complex eigenvalues which may otherwise be 
present for the machine. The real negative roots are 
consistent with the normal operating mode of the 
induction motor, illustrated in the torque-speed (T-w), 
curve of Fig. 4. Generally, the motor operates on the 
negative slope portion of the curve, at the intersection 
point of its mechanical load and the T-w curves. The 
negativity of the E2, E3 and E4 roots denotes a stable 
operation in the negative slope region, following a small 
disturbance about the operating point. 

Figure 3 forms the basis of comparison with the 
parametric variations, described next. 

Parametric Variations 

The effect on dynamic stability, due to 50 to 150 
percent of nominal value variations in LA resistance Ra, 
is minimal while that of La is moderate. This finding is 
similar to that reported in Reference [7]. The impact of 
0 to 150 percent change in the control resistor, Rc, is 
found to be similar to that of Ra. This is evident from 
Fig. 5, and consistent with the circuit topology of Fig. 1. 

Ffffirt of tuning capacitance (£ T 1 - Root locus 
migration caused by 50 to 150 percent change in nomi- 
nal CT is shown in Fig. 6. Here and elsewhere in the 
sensitivity results, the nominal eigenvalue is arrowed, 
where necessary. The trend is similar to that discussed 
in Reference [7]. The impact on the piston (MP), 
displacer (MD) and LA circuit (El) roots is significant. 
However, the induction motor roots E2, E3 and E4 
rem ain stationary, due to their highly damped nature. 
Evidently, nominal value of CT yields a near op tim u m 
placement of the MP root on the jw-axis. 

Fffrrt nf parasitic load (Rp) - The variation in Rp is 
likened to an automatic controller in a practical system. 
However, such incremental changes from 25 to 200 
percent of the nominal value do not account for the 
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n onlineari ties inherent in a physical controller. Nonethe- 
less, the trends portrayed by Rp in Fig.7 illustrate the 
potential effectiveness of a properly designed controller 
in maintaining system stability. 

Other than the roots of the induction motor main 
and backward rotating fields, all remaining roots are 
appreciably affected by the parasitic load. The decay of 
the LA circuit root, El, and the induction motor root, 
E4, becomes more damped, as Rp increases in value. 
The optimum Rp is the nominal value beyond which the 
displacer root MD loses some margin of stability, and 
the piston root MP crosses the jw-axis to become 
dynamically unstable. At the optimum Rp, by virtue of 
the piston root on the jw-axis and all others in the left- 
half plane (LHP), the FPSE/LA subsystem behaves as 
an oscillator. 

Effect of stator short circuit - Typically, system short 
circuit is not a planned event. However, it can happen 
inadvertently. Figure 8 shows the impact of such an 
event on the system stability. Compared to Fig. 3, it is 
evident that the short circuit forces the induction motor 
roots onto the origin. Also, the stability margins of the 
FPSE/LA roots are reduced, as the roots move closer 
towards the right half plane (RHP). The dominance of 
the control resistor R cl in the remaining circuit follow- 
ing the short circuit, assures a near oscillator and stable 
mode of operation. The value of R c j should be chosen 
to ensure system stability at fault condition, while 
minimizing the dissipative loss, due to R c i, during 
normal operation of the power system. 

During the fault condition, the parasitic load Rp is 
effectively shorted and, as such, has no effect on the 
system roots behavior. 

Effect of slip (S) - Variations in motor slip are 
directly related to changes in the induction motor rotor 
speed. A higher slip implies a lower rotor speed, and 
vice versa. An operating induction motor will slow 
down, if its mechanical load, such as a fan, is increased. 
The net result will be a higher slip value. The converse 
is true. The operating slip is determined by the point of 
intersection of the torque-speed and the mechanical load 
curves, namely, point A or B in Fig. 4. 

Figure 9 depicts the system root migration due to 
increasing the motor slip from 25 to 300 percent of the 
nominal value (0.05). Slip has negligible effect on the 
dynamics of the FPSE/LA subsystem. This is partly 
attributable to the relatively small size induction motor- 
fan load (less than liWe), compared to the 12k We 
FPSE/LA subsystem. Also, Fig. 9 shows that the motor 
stator main field root, E2, is independent of slip. The 
effect on the rotor forward field root, E4, is relatively 
more pronounced than that on the backward field root, 
E3. Higher slips reduce the stability margin of E4, and 


slightly increase that of E3. These observations can be 
verified from e xamining the forward and backward field 
circuit impedances, Zf and Zb, respectively, in Fig. 1. 

Effect of no load operation - An excited induction 
motor without mechanical load will run at a slip value 
based on the design requirements of the motor. Fig. 10 
shows the system roots at sudden disconnection of the 
fan load from the motor. Comparison with Fig. 3 of the 
base system shows that the rotor backward field root, 
E3, collapses onto the origin, thereby losing its dynam- 
ics. However, all the remaining roots remain virtually 
unaffected. Thus, for the relatively small induction 
motor /fan load, loss of the mechanical load is not 
expected to substantially affect the dynamic stability of 
the rest of the system. 


Conclusions 

This paper describes the control and dynamic 
stability of FPSE/LA-induction motor /fan system. The 
analysis applies modern control theory to a combined 
dynamic model of the FPSE/LA and circuit and electro- 
dynamic model of the motor/fan load to compute 
eigenvalues of the system. The controllability of the 
system is assessed from the effects of parametric varia- 
tions on the eigenvalues. The key results are summa- 
rized below. 

A system fault, such as a load-side short circuit, can 
potentially force the roots of the dynamic load on to the 
origin and reduce stability margin for the FPSE/LA 
subsystem. The presence of a control resistor on the 
FPSE/LA circuit side promotes a desired oscillatory 
behavior of the engine. A sudden increase in the 
mechanical load beyond the rating of the induction 
motor can reduce the stability of the motor. Depending 
on the relative sizes of the FPSE/LA and the connected 
load, loss of the mechanical load during steady-state 
operation may not substantially affect the system perfor- 
mance. 

The overall results indicate that, with proper design 
and operating conditions and controls, a power system 
comprising a FPSE/LA driving a dynamic load can 
remain dynamically stable. 
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| Subsystem 


Description 


«D 

1.701 (kg) 

Displacer mass 


■v 

9.968 (kg) 
164.63 x IIP (m^) 

Piston mass 


£ 

Piston area 


A> 

102.61 X 10? (tn 2 ) 

Displacer area 


aJ 

2.887 x 10? (m*) 

Piston rod area 


c 

8.1197 x 10 5 (N/n) 

Displacer total stiffness coefficient 



54.95 (N.S/a) 

Displacer total damping coefficient 



12.205 x 10 s (N/n) 

Piston total stiffness coefficient 



60.85 (N.S/o) 

Piston total damping coefficient 



0 (N/m) 

Displacer-piston coupling stiffness co- 
efficient 


K>e 

0 (N/m) 

Piston-displacer coupling stiffness co- 
efficient 



0 (N.S/a) 

Displacer-piston coupling damping coef- 
ficient 


Cpr 

0 (N.S/a) 

Piston-displacer coupling damping coef- 
ficient 

FP Stirling 


186.562 x 10* {N/w 3 ) 

Compression space pressure change with 

Engine 

dXp 

ax,, 

a£ 

-24.98 x 10 6 (H/ra 3 ) 

piston position 

Compression space pressure change with 
displacer position 


0 (N.S/m 3 ) 


Compression space pressure change with 


Up 

« 

0 (N.S/n 5 ) 

piston velocity 

Compression space pressure change with 


0 {N/m 5 ) 

displacer velocity 


dPnt 

Heat exchanger pressure change with pis- 


axp 

0 (N/m 5 ) 

ton position 


2£a 

Heat exchanger pressure change with 



0 (N.S/m 3 ) 

displacer position 


* 

Heat exchanger pressure change with 
displacer velocity 

0 (N.S/m 3 ) 



Heat exchanger pressure change with pis- 
ton velocity 



45.5 (V-S/m) 

Magnetic flux linkage change with piston 
position 

Linear 

R, 

0.082 (Q) 

LA winding resistance 

Alternator 

L, 

1.95 (mH) 

LA leakage Inductance 


C 

6.50 (nfl) 

LA magnetizing Inductance 


Ct 

296.3 (pF) 

Tuning capacitance 


R, 

1.6655 (0) 

Parasitic Load 


Rc 

0.082 (Q) 

Control Resistor 

Single- 

r. 

1.59 (0) 

Stator main winding resistance 

Phase 

Induction 

«-u 

5.2 (*H) 

Stator main winding leakage Inductance 

Motor 

L. 

89.9 (mH) 

Magnetizing Inductance 



6.9 (mH) 

Rotor leakage Inductance (referred to 
stator) 


r'r 

1.869 (Q) 

Rotor resistance (referred to stator) 


S 

0.05 

Operating slip 


Parameter 



Definition 


Ms + MnMr^Mn * ^ 

W<Lm + L,'r) 

L m (0.5L m + ^) + L ls (L m + L| r ) 

Hs * L m L lr * ^ir) 

+ '-tna) 

i; ( L a + Wna) 

Rp/DENOM 

(r $ + R p )/DENOM 
rfTNUM/(2’S*DENOM) 

S*P6/(2-S) 

TNUMT4 

TNUMT5 

TNUM*P6 

-rf*SUM/((2-S) (L m + L^) 2 *DENOM) 

(2-S) *P1 1/S 
S*PlO/(2-S) 









REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information ’Operations and Reports, 1215 Jefferson 
Davis Highway. Suite 1204, Arlington. VA 22202-4302, and to the Office of Management and Budget. Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1. AGENCY USE ONLY (Leave blank) 


2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

February 1994 Technical Memorandum 


4. TITLE AND SUBTITLE 


5. FUNDING NUMBERS 


Controllability of Free-Piston Stirling Engine/Linear Alternator 
Driving a Dynamic Load 


6. AUTHOR(S) 


WU-590-13-11 


M. David Kankam and Jeffrey S. Rauch 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESSEES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-8488 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA TM- 106497 


11. SUPPLEMENTARY NOTES 

Prepared for the 28th Intersociety Energy Conversion Engineering Conference cosponsored by ANS, SAE, ACS, IEEE, ASME, AIAA, and AIChE, 
Atlanta, Georgia, August 8-13, 1993. M. David Kankam, NASA Lewis Research Center and Jeffrey S. Rauch, Sverdrup Technology, Inc., Lewis 
Research Center Group, 2001 Aerospace Parkway, Brook Park, Ohio 44142 (work funded by NASA Contract NAS3-25266). Responsible person, 
M. David Kankam, organization code 5430, (216) 433-6143. 


12a. DISTRIBUTIDN/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 

Unclassified - Unlimited 
Subject Category 66 


13. ABSTRACT (Maximum 200 words) 

This paper presents the dynamic behavior of a Free-Piston Stirling Engine/linear alternator (FPSE/LA) driving a 
single-phase fractional horse-power induction motor. The controllability and dynamic stability of the system are 
discussed by means of sensitivity effects of variations in system parameters, engine controller, operating conditions 
and mechanical loading on the induction motor. The approach used expands on a combined mechanical and thermo- 
dynamic formulation employed in a previous paper [Ref. 7]. The application of state-space technique and frequency 
domain analysis enhances understanding of the dynamic interactions. Engine-alternator parametric sensitivity studies, 
similar to those of the previous paper, are summarized. Detailed discussions are provided for parametric variations 
which relate to the engine controller and system operating conditions. The results suggest that the controllability of a 
FPSE-based power system is enhanced by proper operating conditions and built-in controls. 


14. SUBJECT TERMS 

Free-piston Stirling engine; Dynamic load; Dynamic stability; 
Controllability; State-space technique 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 

NSN 7540-01-280-5500 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


15. NUMBER OF PAGES 

9 


16. PRICE CODE 

A02 


I 20. LIMITATION OF ABSTRACT 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 















